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We have obtained an analytic expression for the ring diagrams contribution to the correlation 
energy of a two dimensional electron liquid as a function of the uniform fractional spin polarization. 
Our results can be used to improve on the interpolation formulas which represent the basic ingredient 
for the construction of modern spin-density functionals in two dimensions. 
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The electron liquid, an ideal model in which electrons 
interact via a Coulomb potential in the presence of a 
uniform rigid neutralizing background, represents a fun- 
damental paradigm for the understanding of condensed 
matter systems^ Although for homogeneous states the 
ground state properties are completely determined by the 
number density n, it can be useful and often necessary 
to study the properties of the system in the presence of 
a given uniform spin polarization. Accordingly, for the 
case under consideration of a two dimensional system, 
alongside the traditio nally em ployed dimensionless den- 
sity parameter r s = y/ ima 2 B we will consider the effects 
of polarization as parameterized by p — , where S z is 
the uniform spin polarization density. 

While in the general uniform case the total energy of 
the system can be accurately obtained only by numerical 
means and is available in the form of interpolation for- 
mulas of Monte Carlo resultsji^ an exact analytic treat- 
ment is possible in the high density limit in the form of a 
perturbative expansion in r s which in the case is vanish- 
ingly small^ In this limit the first-order correction to the 
nonintcracting energy simply coincides with the famil- 
iar exchange energy, which has a simple analytic formula 
proportional to rj 1 . The remaining correction is referred 
to as the correlation energy^ 

The high-density expansion of the correlation energy 
was studied in the classic Ref. 0- The final expression for 
the paramagnetic case (in Ry units) reads: 

£ c (r s ,p = 0) =-0.385-^p(10-37r)r s lnr s +... , (1) 

while the p = 1 result can be also obtained, by making 
use of a simple transformation^ 

For generic values of the polarization, the coefficients 
of the expansion (fT]) become functions of p. 

The constant term in |T]) is obtained from second-order 
perturbation theory and is the sum of two distinct con- 
tributions. The first is the second-order exchange en- 
ergy, which is independent of p and can be calculated 
analytically.— The second one stems from the second- 
order direct energy term, commonly referred to as the 
first ring diagram. The latter has been recently evalu- 
ated numerically for generic values of p in 

Rcf.ffl These 

contributions are accurately represented by the interpo- 



lation formula of Ref. 0. 

The next perturbative terms are in general divergent 
and a finite result is obtained upon exact summation of 
the infinite series of the most diverging contributions, the 
remaining ring diagrams. This leads to the sub-leading 
r s In r s term. 

This elegant method was originally developed for the 
corresponding unpolarized three dimensional case^ a 
problem in which the ring diagrams sum up to give the 
leading In r s contribution to the correlation energy £ For 
the three dimensional case the exact polarization depen- 
dence was also determined. 9 

The generic dependence of the r s In r s term on p is 
obtained in this work. The complete formula for the 
presently known leading terms of £ c {r s ,p) is provided 
for reference in the Appendix. 

Following Ref. 0, the value of the generic diverging ring 
diagram of order n is obtained from the expression: 

An), v (-1)" [ + °° , [°° 2h ( Q q (u)r s \ n 

(2) 

where n is a positive integer larger than 1. 
The explicit expression for Q q {u) is given by: 

<M«) = JI / 7 (f + ^ M*)(i-M*0), 

(3) 

where k' = ^J{k x + q) 2 + k 2 . The polarization depen- 
dence of Q q (u) is implicitly determined by the occupation 
functions: 

n a {k)=e{k a -k) . (4) 

For each of the two spin orientations, the Fermi wave 
vector is obtained from the relation: 



fc T(i) = V l± P > ( 5 ) 

where for convenience we have rescaled the wave vectors 
by kF — y/2nn. 

The expression corresponding to the generic ring dia- 
gram of Eq. ([21) has a (formal) dependence of r™~ 2 , but 
in reality for n > 3 it involves a diverging integral. The 
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sum to infinite order is 
1 



/fl(p) = 



du 



q 2 dq 
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which is instead a converging integral whose leading 
contribution in the r s — > limit is again proportional 
to r s In r s and is determined by the small-g integration 
region. 10 Accordingly it is sufficient to employ here the 
limiting value Q q =o( u ): a s obtained from ([3]). The result 
can be written as: 



Qo(u) = 2n[R(u/k ] ) +R(u/k{)] , 
where the function R(u) is defined as: 

R(u) = 1 1 . 

y/1 + l/v 2 



(7) 



(8) 



To this point our discussion follows Ref. 0, where the 
implicit expression ([6]) as well as Eq. ((7|) were originally 
provided for a generic value of p. There however, the 
explicit evaluation was only done forp = 0. In the general 
case, using ([7]) and performing in ([5]) the wave vector 
integration (up to an arbitrary upper limit), one extracts 
the leading contribution in r s : 



f R (p) 
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] {u)] 3 du , 



(9) 



which gives the standard result |T]) by making use of 
Qo(u) = 4itR(u), i.e. 
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We find that the integral in ([9]) can be performed ex- 
actly in the general case. The result can be expressed in 
terms of the function 



F(x,y)= / [R{u/x)] 2 R(u/y)du 

J — OO 

which has the explicit expression: 



(ii) 



y 2 , 2x 2 arccos - 

F(x, y) = 4(a + y)-7rx- 4xE(l - 1L) + f . 

x ^Jx l - y z 

(12) 

Here E(x) is the complete elliptic integral of the sec- 
ond type,— and one should use the identity arccos * = 

\/x 2 -y 2 

arccosh x ^ y > x. At p = one only needs the value 

F(l,i) = 10- 3tt. 

The result for the ring diagrams at finite polarization 
can be compactly and elegantly expressed in terms of 
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FIG. 1: Plot of the scaling factor I R (p) (defined in Eq. fl3)l) 
as a function of the fractional polarization p. Solid line: exact 
expression (|14[) . Dashed line: Eq. (|15p as obtained from the 
interpolation formula proposed in Ref. |2|. 



the corresponding polarization scaling function Ir(j>), de- 
fined as: 

Ir(p) = lim n TW ' ( 13 > 

The final expression is given by: 

If, , F(kt,k\) + F(k i) kt)\ , > 
Jh(p) = g Ut + fc l + 3 T 1 T - j . (14) 

which readily gives the correct value at p = 0. At p = 1, 
using F(\/2, 0) = F(0, \/2) = 0, we obtain the known 
result* I R (1) = ^. 

As it turns out the exact result Ir(p) is not well 
reflected in the most recent interpolation formulas of 
Monte-Carlo calculations provided in the literature. In 
particular, from the correlation energy formula of Ref. 2, 
denoted here as £jr (r s ,p), the following limit is ob- 
tained: 

e lc (r s ,p)~£* IC (0,p) 



I^(p)= lim 



M0) 



1 - 0.3932p 2 - 0.4297p 4 



(15) 



which is compared in Figure Q] to the exact result of 
Eq. (|14p . The difference is remarkable, even if the specific 
aim of Ref. [2j is to address the polarization dependence of 
the whole correlation energy. Agreement of ijj (p) with 
the exact result is only achieved for p — and p = 1, 
values known from the extant literature. A noticeable 
failure is the behavior near p = 1 , where the polynomial 
(fT5f gives a finite slope while the leading term in the 
exact expression is: 



T ( \ ^ 



14- 3tt 



(16) 



4(10-37r) 

Around p = 0, while being correctly quadratic in p, 
Eq. (fl~5j) displays an incorrect coefficient. The exact co- 
efficient is given by: 

168 - 45tt 



I Rip) 



1 



160(10 -3tt) 



P 



1 - 0.2893 p 2 



(17) 
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The disagreement between I(p) and the Monte Carlo 
based interpolation I M (p) is hardly surprising for the 
latter was obtained by sampling the energy at a num- 
ber of polarization values for each of the values r s — 
1,2,5,10, which are clearly outside of the r s -C 1 per- 
turbative regime. 12 In practice, I(p) refers to the sub- 
leading term in the density expansion, so that it gives 
only small corrections to the total energy. Nevertheless, 
incorporating this exact formula would certainly result 
in an improved empirical expression for the polarization 
dependence of the correlation energy. 

We conclude by noting that a similar, yet more in- 
volved, calculation can be carried out for the high density 
limit of the electron liquid correlation energy in the pres- 
ence of Rashba spin orbit. An analysis of this interesting 
and timely problem is provided elsewhere^ 



APPENDIX 

For ease of reference, we collect here the explicit form 
of all the leading terms contributing to the perturbative 
expansion of the total energy of the two dimensional elec- 
tron liquid at finite polarization. The general formula (in 
Ry units) reads: 

£(r s ,p) = £ K (r s ,p) + £ x (r a ,p) + £ c (r s ,p) . (A.l) 

In this expression £k{t s ,p) represents the noninteracting 
kinetic energy and is given by: 



while £ x (r s , p) represents the exchange energy: 

£l(r ,, p) = _^(i±iM^. (A.3) 

The expansion for the correlation energy takes the form: 
2 v/2 

£ c (r s ,p) = £2(p) 7, — (10 - 37r) Jr(p) r s \nr s + ... , 

where the omitted corrections are of order 0(r s ). In this 
expression Ir(p) is defined by Eqs. (fT4|) . (fl"2f and ([5]). 
The density independent term is in turn given by: 



where ^ 



and 6 



£ 2 (p) = 4 b) +4 r \p), 



£^ b) = 0.2287 , 



4 r) (p) = -0.6137 1 2 {p) , 



(A.5) 



(A.6) 



(A.7) 



p < \ 1+ P 2 



where the scaling function l2{p) is given by: 

(l+p)ln(l+p) + (l-p)ln(l-p) Sf(p) 

h{p) = l 4hr2 2~' 

(A.8) 

witbM 

Sf(p) ~ 0.0636p 2 - 0.1024/ + 0.0389/ . (A.9) 



(A.2) 
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